Magnetic impurity in the vicinity of a vacancy in bilayer graphene 
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We use quantum Monte Carlo method to study a magnetic impurity located next to a vacancy 
in bilayer graphene with Bernal stacking. Due to the broken symmetry between two sublattices in 
bilayer system, there exist two different types of vacancy induced localized state. We find that the 
magnetic property of the adatom located on the adjacent site of the vacancy depends on whether 
the vacancy belongs to A or B sublattice. In general, local moment is more strongly suppressed if 
the vacancy belongs to the sublattice A when ^ ~ 0. We switch the values of the chemical potential 
and study the basic thermodynamic quantities and the correlation functions between the magnetic 
adatom and the carbon sites. 



I. INTRODUCTION 

Interest in the defects induced localized states in 
monolayer and multilayer graphene has been greatly 
expanding^^^ driven by the enthusiasm to seek for ferro- 
magnetism in graphene^^^. Experimentally, one can in- 
troduce vacancies in carbon-based systems by irradiation 
to modify their properties'^ In monolayer graphene, 
vacancy leads to the formation of zero-energy quasi- 
localized states^ such that the local density of states 
(LDOS) of the neighboring sites would be greatly mod- 
ified. Consequently, the behavior of a magnetic impu- 
rity in the vicinity of a vacancy shall be different either 
from that in a pure graphene or from that in a nor- 
mal metal. Study on the consequences of an adjacent 
vacancy-adatom^^ showed that the local moment forma- 
tion is strongly suppressed. However, in pure graphene 
where the density of states(DOS) vanishes at Dirac point, 
the Kondo effect does not occur in general so a well- 
developed local moment will show up in undoped case^^. 

In Bernal stacked bilayer graphene (BLG), due to the 
interlayer hopping energies, the two sublattices are no 
longer equivalent. It has been reported that whether the 
defect is generated on A or B sublattice would greatly 
alters the defect induced magnetism and quasi-localized 
states^ 

In this paper, we study the magnetic properties of 
an Anderson impurity on the adjacent site of a va- 
cancy in Bernal stacked bilayer graphene using the 
quantum Monte Carlo method based on the Hirsch-Fye 
algorithm^^. 

We invoke the Anderson impurity model, and the total 
Hamiltonian can be written as 

H = Ho + Hi+H2+ H^ac (1) 

Hq, Hi, H2 are the tight-binding Hamiltonian of BLG, 
the impurity Hamiltonian and the hybridization between 
impurity and the carbon electrons, respectively. They 



can be written as 

Ho = {-t ^ aL^"u - ^1 XI "L«2i) + H.c, 

<i,j>m j 

Hi = ll-^ialmO-rm + bl^.fim^), (2) 
in I 

Hi = {Ed - li)dU + + H.c), 

where amiaibmia) annihilates an electron with spin a at 
the site Rmia {Rmib) on sublattice A(B) of graphene's 
hexagonal structure in the m-th layer, we use t k, 
2.8eF^^ as the energy unit. For the inter-layer hopping 
energies, according to the Slonczewski-Weiss-McClure 
parametrization-'^^'^', ti w Q.AeV . annihilates an elec- 
tron with spin a at the impurity orbit. If the magnetic 
adatom is added on the top of sublattice A, Cmda = oLmdai 
otherwise Cmda = ^mdo-, where d denotes the location of 
the impurity adatom. Without loss of generality, we as- 
sume the vacancy is located on layer 1. If the vacancy is 
adsorbed on sublattice A, the Hamiltonian Hvac induced 
by the vacancy can be written as 

Hvac ^'t'^{a\obij + H.c.) + ti{a\^a2a + H.c.) 

(3) 

We would not consider about the effect of reconstruction 
of the remaining structure in the presence of a vacancy, 
since the properties of localized states are rather insensi- 
tive to the deformation-'^^. 



II. RESULTS 

We mainly test two different groups, namely, vacancy 
belongs to sublattices A and B, and the magnetic adatom 
is located on the adjacent site of the vacancy. In all the 
results we show below, we apply the minimal model and 
use ti = Q.2t and ignore other interlayer hopping energies 
since they are much smaller than ti and have minor effect 
on the behavior of magnetic adatom. 
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FIG. 1: (Color online). The various thermodynamic quanti- 
ties for the two cases that vacancy embedded in sublattices 
A and B. The impurity is located on the nearest site of the 
vacancy on layer 1. (a) Ud as a function of ^. (b) Udfridi as a 
function of /i. (c) m| as a function of /i. (d) x as a function 
of fi. In all the plates we use V = l.Ot, U — 0.8t, Sd = ~U /2 
and /3 = 1/r = AOt'^ . 



Shown in Fig. 1 are the basic thermodynamic quanti- 
ties of the impurity on the nearest neighbor of the va- 
cancy for both the cases that the vacancy belongs to 
sublattices A and B. Shown in Fig. 1 (a) are the val- 
ues of charge Ud, and we can see that the A vacancy case 
shows slightly larger occupancy than B vacancy case near 
the zero-energy. The values of double occupancy Ud^Udi 
shown in Fig. 1(b) has the same order as that of Ud- 
We can see that as we lower the values of chemical po- 
tential from zero-energy, the values of Ud and nd-\nd\. are 
gradually decreased. However, the local moment squared 
shown in Fig. 1(c) is first increased and then decreased 
for the A vacancy case, while mostly unchanged for B 
vacancy case around the zero-energy. According to the 
fact that nda can either be zero or one, we have 



md = nd- 2nd^nd\,. 



(4) 



The suppression of the local moment formation is mainly 
caused by the larger values of double occupancy Ud^Udi 
compared with the pure BLG case. We can also find that 
the impurity magnet moment for the B vacancy case is 
slightly larger than that of A vacancy case. As the chem- 
ical potential is lowered, they cross over each other and 
finally the difference between these two values vanishes. 
The difference is more obvious in the values of spin sus- 
ceptibility as shown in Fig. 1(d). In principle, spin sus- 
ceptibility depends not only on the local moment itself 
but also on the spin correlation with the conduction band 
electrons. The A vacancy leads to sharper resonance on 
the LDOS of the nearest neighbor around E = 0, which 
will cause stronger screening of the local moment. From 
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FIG. 2: (Color online). (a)The local moment and (b)the sus- 
ceptibility for the two cases that vacancy on sublattice A and 
on sublattice B as a function of fi. The impurity is located 
on the nearest site of the vacancy. In all the plates we use 
V = 0.5t, £d = -(7/2 and ^ = l/T = 40t"\ 




FIG. 3: (Color online). The spin susceptibility as a function 
of T"^ for the case that (a) vacancy on sublattice A and (b) 
on sublattice B. The impurity is located on the nearest site 
of the vacancy. In all the plates we use V = l.Ot, Ed = —(7/2. 



Fig. 1(d), we can also find that the spin susceptibility is 
first increased and then decreased as we lower the chem- 
ical potential. Presented in Fig. 2 are the values of 
and X for y = 0.5t, /? = 40t^^. The values of the local 
moment shown in Fig. 2(a) and the susceptibility shown 
in Fig. 2(b) arc much larger than those for V = l.Ot case 
as shown in Fig. 1. However, the general behavior of the 
local moment and the spin susceptibility with respect to 
/i is the same as the case oiV = l.Ot that for both A and 
B vacancy, the and x ^-re first enhanced and then 
suppressed as the chemical potential is lowered. Fig. 3 
shows the spin susceptibility of the adatom as a function 
of inverse temperature for V = l.Ot and £d = —U /2. For 
both cases, we see that at low temperature, the suscep- 
tibility is temperature-independent at = 0. As /i is 
lowered, we can see that the values of x tend to regain 
the features of Curie-like dependence. The reason is that 
as ^ is moved away from the zero-energy, the LDOS of 
the carbon sites near the vacancy drops significantly and 
as a result, the screening of the local moment is sup- 
pressed. Hence, a well-developed local moment shows up 
on the impurity site. However, if we go on lowering /i, 
a temperature independent behavior would be regained. 
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the charge-charge correlation and spin-spin correlation 
would decrease in values. We can also find that both 
the spin-spin and charge-charge correlations between the 
magnetic impurity and the carbon site show larger ampli- 
tude for the A vacancy case. This is consistent with the 
results of the basic thermodynamic quantities presented 
in Fig. 1. 



III. SUMMARY 



FIG. 4: (Color online), (a) The charge-charge correlation Co 
and (b) the spin-spin correlation with respect to ^. For 
both cases the magnetic adatom is attached to the nearest 
site of the vacancy. We use V = l.Ot, U = 0.8t, Ed ~ —U/2 
and /3 = 1/T = 40i"^ 

As we further lower the chemical potential, charges as 
well as the spins move out from the impurity and carbon 
sites, leading to a decrease in the spin susceptibility. 

In Fig. 4, we plot the behavior of on-site correla- 
tion functions Co and So with moving below the zero- 
energy. Here Co and 5'o is the correlation between the 
magnetic impurity and the carbon site where the impu- 
rity sits. We use V = l.Ot, U = 0.8t, Sd = -U/2 and 
/3 — 1/T = 4:0t~^. We can sec that both the spin-spin 
and charge-charge correlations decrease in values as the 
chemical potential is lowered, as is shown in Fig. 4. Both 
the correlations are the strongest at half-filling, and as 
the chemical potential is lowered, the correlations are 
reduced in amplitudes. As we have analyzed, vacancy 
induces zero-energy localized states on the surrounding 
sites and hence, at /.i = 0, the impurity spin is strongly 
screened by the conduction electrons in the host material. 
As the chemical potential is lowered, both the impurity 
site and the carbon sites would lose the charge such that 



In summary, we studied an on-site magnetic impurity 
placed in the vicinity of a vacancy in the bilayer graphene 
with Bernal stacking. We mainly test two different cases, 
namely, vacancy located on A and B sublattices. We find 
that if the magnetic adatom is located on the adjacent 
site of the vacancy, the local moment is strongly sup- 
pressed due to the zero-energy mode which is induced by 
the vacancy. The local moment for the A vacancy case 
are slightly smaller than that in B vacancy case at /i ~ 0, 
due to the broken symmetry caused by the Bernal stack- 
ing. As /i is lowered, both the spin and charge correla- 
tions between the adatom and carbon sites are weakened. 
This behavior is completely different from that in pure 
BLG, where both of the correlations arc enhanced as the 
chemical potential is lowered. 
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